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Abstract 



We show that the algebra of discretized spatial diffeomorphism constraints in Hamiltonian 
lattice quantum gravity closes without anomalies in the limit of small lattice spacing. The 
result holds for arbitrary factor-ordering and for a variety of different discretizations of the 
continuum constraints, and thus generalizes an earlier calculation by Renteln. 



Motivation 



For the most part, our understanding of quantum gravity in four dimensions is not 
sufficiently developed to provide unambiguous tests for the consistency and appropriateness 
of any candidate theory. In the case of usual gauge field theories, there are well-known criteria 
to ensure their consistency at the quantum level, for example, the validity of the Slavnov- 
Taylor identities. However, the Fock space methods on which they are based are only of 
limited use for gravitational theories in more than two space-time dimensions. 

Attempts to construct non-perturbative path integral formulations of 4-d gravity have 
so far left the question of their continuation to Lorentzian signature unanswered. Even dis- 
regarding this problem, it is difficult to get a mathematical handle on the structure of the 
path integral measure. Canonical approaches circumvent the signature problem and allow 
one to pose (part of) the quantum consistency problem in a seemingly clear-cut form: is there 
an anomaly in the Dirac algebra of the Hamiltonian and diffeomorphism constraints or isn't 
there? In the presence of anomalous terms, the quantum Dirac conditions overconstrain the 
space of physical states, which is physically unacceptable. One difficulty of the canonical for- 
mulation lies in finding appropriate regularized versions of the constraint operators, without 
which the question is known to be ill-posed [1]. A potential drawback of the Hamiltonian 
approaches is the fact that efficient computational techniques for a numerical study of their 
quantum properties have yet to be developed. 

Schemes that proceed by a direct discretization of space-time or of the space of all 
field configurations in order to achieve a regularization tend to break the diffeomorphism 
invariance present in the continuum theory. This however is not a fundamental objection, if 
one considers the smooth structure of the classical theory only as a semi-classical property, 
which does not continue to hold down to smallest length scales. Still it leaves one with the 
technical inconvenience of having to work with an "approximately diffeomorphism-invariant" 
regularized theory. 

Considerable efforts have gone into trying to construct a canonical quantization of 3-1-1 
gravity in terms of connection variables, without resorting to a discretization, and with the 
full spatial diffeomorphism group still acting on a suitably defined space of quantum states 
[2-5]. However, this approach seems to be running into some difficulties, among which are a 
"non-interacting" property of typical realizations of the quantum Hamiltonian, which leads, 
for example, to the existence of unexpected local quantum observablcs [6,7], vanishing com- 
mutation relations between two such Hamiltonians [8] (even in situations where the diffeo- 
morphism group acts non-trivially) , and anomalous commutation relations between geometric 
operators that classically commute [9,10]. 
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This provides an additional incentive for studying a truly discretized version of this 
formulation, which does not share any of these features. On the other hand, like most 
discretized models, it does not carry any obvious, non-trivial representation of the classical 
invariancc group of general relativity. The diffeomorphism group is only recovered in the 
limit as the regulator is taken to zero. 

In the present piece of work, we will generalize a result obtained earlier by Renteln [11] in 
the framework of a lattice discretization of 3+1 gravity in a connection formulation [11-17]. 
Various aspects of this version of lattice quantum gravity have changed since these earlier 
investigations, as a result of new developments both in the continuum loop representations of 
quantum gravity mentioned above, and in lattice gravity proper. The kinematic resemblance 
with Hamiltonian lattice gauge theory has become even closer since the theory was rephrased 
in terms of real su(2)-valued connections [18], instead of the original sl{2, C)-valued Ashtekar 
potentials [19], in order to avoid complications related to the reality conditions that appear 
in the complex case. As a by-product, there is now a well-defined scalar product and a 
complete set of square-integrable states for each finite-size lattice. Although this may not 
give rise to a scalar product on the sector of physical states (satisfying all of the constraints), 
it is nevertheless remarkable that it should exist at all, given that inner products are usually 
hard to come by in quantum gravity. 

We will show that independent of the factor- ordering and for a variety of operator sym- 
metrizations, there are no anomalous terms in the commutator algebra of the lattice analogues 
of the spatial diffeomorphism generators in the limit of vanishing lattice spacing. The compu- 
tations are already substantial for this subalgebra of the entire quantum constraint algebra. 
We will also spell out some details of the calculations, that may be of interest in attempts 
to compute the commutator in alternative regularization schemes. This result shows that 
no inconsistencies arise in the lattice discretization at this level, and paves the way for the 
calculation of commutators involving also the Hamiltonian constraint. 



1 Introduction 

A prominent feature of the lattice discretization we are considering is its resemblance with 

canonical lattice gauge theory. Like Hamiltonian gauge field theory, 3-|-l continuum gravity 
can be formulated in terms of Yang- Mills conjugate variable pairs ( A, E) of "Ashtekar-type" , 
with Poisson brackets {A\{x),E^j{y)} = y). The variable is a spatial stt(2,IR)- 

valued gauge potential. 
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Up to terms proportional to the Gauss law constraint, the Hamiltonian constraint in this 
formulation is the Barbero Hamiltonian [18] (rescaled by a factor (det£')~2 to make it a 
density of weight one), plus - for the sake of generality - a cosmological constant term. 



n{x) =-^e^^^EtE]F^, - ^{detE)-'/'r,a,asaMsbAEt'E^-E'^E'^- 
- 2E^E:^El^E^')El^E'^iVa,En{V,,E'^') + ^Vd^ = 0. 

We have chosen the dimensions of the basic variables and constants to be [A] = L~^, 
[E] = , [G] = L^, [A] = L~'^ , where G is Newton's constant. Also recall that E transforms 
as a one-density under spatial diffeomorphisms. In (1.1), TL has already been brought into 
the form of a polynomial modulo powers of \/deLE, in order to make its discretization and 
quantizion straightforward. The spatial diffeomorphism constraints are given by 



V,{x) = Et{x)FUx) = 0. (1.2) 

Smearing out these four constraints with the lapse and shift functions N{x) and iV"(x), one 
arrives at the expressions 

H[N] ■= f Sx N{x)H{x) 

^ r (1-3) 

y[iV"] := / d^x N^{x)Vb{x), 



which satisfy the usual Dirac Poisson algebra. In particular, for the spatial diffeomorphism 
generators one derives 



{V[N''], V[M'']] = j d^x (N'^dbM'' - M^dbN'')Fi^El - J d^x F^,,N''M\VcE^) 

which coincides with the diffeomorphism Lie algebra on the subspace of phase space given 
by the vanishing of the integrated Gauss law constraints G[A*] = f d'^x A^{\7 aEf). One may 
also redefine the generators y[A^"] by adding a suitable term proportional to the Gauss law 
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constraint in order to get rid of the extra phase-space dependent term on the right-hand side 
of (1.4) (see, for example, [20]). 



2 Discretization 

Before setting up the discretization of operators relevant for lattice gravity, we start 
with a brief summary of the basic ingredients of Hamiltonian lattice gauge theory [21]. For 
computational simplicity, we take the lattice A to be a cubic AT^-lattice with periodic bound- 
ary conditions. The cubic symmetry is convenient but not necessary, as we know from the 
continuum loop quantizations which use imbedded lattices. 

The basic quantum operators associated with each lattice link I are a group- valued SU (2)- 
link holonomy IJ (represented by multiplication by U), together with its inverse U~^, and a 
pair of canonical momentum operators pf and , where i is an adjoint index. The operator 
pf{n,a) is based at the vertex n, and is associated with the link I oriented in the positive 
a-direction. By contrast, p^{n -|- la, a) is based at the vertex displaced by one lattice unit in 
the a-direction, and associated with the inverse link l~^{a) = l(— a). In mathematical terms, 
the momenta p"*" and p~ correspond to the left- and right-invariant vector fields on the group 
manifold associated with a given link. The wave functions are elements of XiL'^{SU{2),dg), 
with the product taken over all links, and the canonical Haar measure dg on each copy of the 
group SU(2). The basic commutators are 



[t/A^(n,a),Z7c^(m,6)] =0, 

\pf{n,a),UA'^{m,b)] = -^SnmS^i,TiA^UB'^ {n,a), 

\p~{n,d),UA^{'m,b)] = -^Sn,m+iS^f,UA^ {n,d)TiB'^ , 

[Pj^(^,«),Pj'("i,^)] = ^iSnmS^i,eijkPk{n,a), 
\pf{n,d),p~{m,b)] = 0, 



where eijk are the structure constants of SU{2). (There are analogous Poisson bracket re- 
lations for the corresponding classical lattice variables.) The commutation relations for the 
inverse holonomy operators can be easily deduced from (2.1). The normalization for the 
S'?7(2)-generators is such that [Ti,Tj] = 2ejjfeTfe, Tr TjTj = —2Sij and = Al^Ti/2. 

In order to relate discrete lattice expressions with their continuum counterparts, one uses 
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power series expansions in the so-called lattice spacing a, which is an unphysical parameter 
with dimension of length. For the basic classical lattice variables, these are 

pf(b) = El ± G'^ VbE^ + O(a^). 

We will assume that similar expansions continue to be valid in the quantum theory. Note 
that Newton's constant G appears in (2.2) since the dimensions of the basic gravitational 
variables A and E differ from those of the corresponding Yang-Mills phase space variables. 

Using the expansions (2.2), one obtains unambiguous continuum limits of composite 
classical lattice expressions by extracting the coefficient of the lowest-order term in the a- 
expansion. The converse is not true: there is no unique lattice discretization of a continuum 
expression, since one may always add to the lattice version terms of higher order in a which 
do not contribute in the continuum limit. We will consider in the following different sym- 
metrizations for the lattice diffeomorphism generators. 



3 Further preparations 

As a first step in the commutator calculation of two lattice diffeomorphisms, one needs to 
expand the classical link holonomies in a neighbourhood of the lattice vertex no at which the 
local commutator will be computed. For simplicity we choose a local lattice coordinate system 
such that the vertex no coincides with its origin, that is, no = (0, 0, 0). There are six nearest 
neighbours, which are one lattice unit away, for example, the vertices n = (±1, 0, 0) in positive 
and negative 1-direction. Because of the non-locality of the diffeomorphism generators, all 
relevant link holonomies lie on a cube of edge length 2 about the origin. 

Using the defining formula for the holonomy along a parametrized path 7 between two 
points xi and X2, with 7(0) = xi and 7(1) = X2, 



00 „i „i „i 
U{xi,X2) = t + Y, dsi dS2... ds„ A(7(si))^(7(s2))...^(7(sn)), (3.1) 
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(2.2) 



one expands this into a power series in the lattice spacing a for the special choice Xi = (0, 0, 0), 
X2 = (O) 0, 0), say. We will need terms up to power for our computation. For the holonomy 
in positive 6-direction based at (0,0,0), one obtains thus 



U{{0, 0, 0), b) = l + aA, + -{d,A, + Al) + -(^^^^ + id,A,)A, + 2A,id,A,) + Al) + 0{a^), 

(3.2) 

where we now have set G = 1. (Note that at order this differs from the expression given 
in [11].) The other positively oriented link holonomies are computed similarly, as Taylor 
expansions with respect to the local cartesian coordinate system spanned by the lattice. From 
those, the inverse holonomies (associated with the opposite link orientation) are computed 
order by order using UU~^ = 1. Expansions for the link momenta are obtained in an 
analogous manner, using (2.2). 

Two useful checks for possible errors are given by calculating holonomies Ua around 
plaquettes, i.e. shortest closed lattice paths of edge length four. Denoting by C/n^^ the 
plaquette that starts with a link in positive 1-direction, followed by links in directions 2,-1 
and —2, one has the expansions 

Tr Ua,, = 2 + 0{a*), TV (Ua^^n) = a^Fl^ + ©(a^). (3.3) 

A natural discretization for the smeared spatial diffeomorphism generators, (second line 
of (1.3)), is schematically given by 

^iVi-"(n,a)Tr(?7n,,r,)p,(n,6) =: ^ V^-"[iV, n). (3.4) 

n n 

where we have not yet specified the details of the implicit summation over the indices a 
and h. We will consider four distinct ways of symmetrizing the local lattice expressions 
A'"^'^**(n, a) Tr (f/n^(jrj)pj(n, 6). Let us introduce the notation N^{n,d) for the lattice shift 
functions at vertex n, associated with the link in positive and negative a-direction respectively. 
This notation conforms with the one for the momenta p^(n, a). In order for (3.4) to have 
the correct continuum limit as a ^ 0, we require that 
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N+{n, b) ^ -N\x) + N^{x) + O(a^) 



1 - (3.5) 

N-{n, b) ^ -N\x) + iV^(x) + O(a^). 

We leave the zeroth-order functions N^{x) and N^{x) unspecified; the final result will not 
depend on this choice. For the local continuum expression of the form N'^F2^Ef, say, we will 
consider four different lattice representations: 

(a) no symmetrization: 

iV+(n,2)TV(i7n,3)p+(n,3); 

(b) symmetrization over shift functions: 



-(iV+(n, 2)Tr {Uo,Min, 3) + N- {n, 2)Tr {Uo,,_Min, 3)); 



(c) symmetrization over momenta: 



^(iV+(n, 2)Tr iUa,s)pt{n, 3) + N+{n, 2)Tr {Ua_s,,)p- {n, 3)); 



(d) symmetrization over both shift functions and momenta: 

J(7V+(n, 2)IV (C/n,3)p+(n, 3) + N-{n, 2)Tr ([/□3,_,)p+(n, 3)+ 
iV+(n, 2)TV (t/n_3 Jp-(n, 3) + N-{n, 2)Tr {Ua_,^_,)p- i^^ 3)); 

Note that in each case we sum over terms that are maximally localized on the lattice, 
i.e. consist only of holonomies and momenta located on a single lattice plaquette. Choice (d) 
was the one considered in [11]. 



4 The commutator 

We now come to the actual computation of the commutator of two lattice-regularized 
diffeomorphism constraints. As is well known, the algebra of the discretized constraints 
(3.4) does not close, not even at the level of the classical Poisson algebra, because the result 
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is a sum of terms, where each term may extend over up to two plaquettes. That is, the 
Poisson commutator is not a Unear combination of terms of the form (3.4). However, if the 
discretization is consistent, the commutator should yield the continuum result to lowest order 
in the lattice spacing a as a — > 0. 

We will perform the corresponding quantum computation, and for definiteness will choose 
initially the factor-ordering for the Tr (Ua^^Ti)pi{n, 6)-terms with all momenta p to the right. 

To obtain the contribution at a given lattice vertex uq to the commutator of two smeared 
lattice diffeomorphisms V[N, n),Y^^ V[M, m)], one has to collect all terms where a mo- 
mentum based at no acts on either a holonomy or another momentum based on the same link. 
One quickly realizes that there are not only contributions from terms in [V[N, no), V[M, no)], 
but also contributions from vertices close-by, either from nearest neighbours (e.g. no±(l, 0, 0) 
- there arc 6 such vertices) or from vertices across a plaquette diagonal (e.g. no ± (1, 1, 0) - 
there arc 12 such vertices). We will see below that these non-local contributions are indeed 
necessary for obtaining the correct result. 

Without loss of generality, we will consider local lattice smearing functions N{no,a) = 
(7V^,0, 0) and M(no,a) = (0, M^,0). For the fully symmetrized version of the local con- 
straints V[N,n), there are 36 non- vanishing commutators between terms based at no and 72 
non-vanishing commutators involving also terms based at neighbouring vertices. The commu- 
tators are evaluated using the basic commutation relations (2.1) and various epsilon-function 
identities. The most time-consuming task is the expansion and subsequent simplification of 
the resulting expressions in powers of a. It turns out that the holonomies have to be expanded 
to third order, i.e. up to terms cubic in the local connection form A(x), and there are rather 
a lot of terms at this order. 

In the appendix we have collected the separate contributions from both the local com- 
mutators and those involving neighbouring vertices, for all four types of symmetrization. It 
is instructive to see what kind of terms arise in the different parts of the calculation and 
how they cancel. This may be relevant to attempts in the continuum loop representation 
of quantum gravity to reproduce similar quantum commutators. Note also that in order to 
derive (4.1) below from the formulae given in the appendix, some partial integrations had to 
be performed. 

After a lot of algebra one finds that, independent of the symmetrization, the final result 
is given by 
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lim[^V[iV,n),^V[M,m)]U = 

n m 

a^{N\diM^){F^sEf - Fl^E}) + M\d2N^){HiEf - Pt^Ef) 

N^M\ViF^^ + VsFa'i + ^3Fh)Ef) + 0{a^). 

Comparing with equation (1.4), this is the expected answer, without any anomalies, and up 
to a term proportional to the Bianchi identity. It shows that at this level both the classical 
discretization and the quantization of the diffeomorphism constraints are consistent. The 
independence of the symmctrization suggests that there is a chance that the evaluation of 
more complicated commutators involving also the Hamiltonian constraint may already yield 
the right result if done in terms of the unsymmetrized lattice constraints. This is potentially 
important since the calculations become even more complicated. 

Let us now turn to the issue of factor-ordering. We have already proven the absence of 
anomalies with momenta ordered to the right and will now use this result to deduce what 
happens for arbitrary factor-ordering. Let us adopt the notation (Tr C/ri)|»i for a typical term 
of the quantized diffeomorphism constraint, where U denotes some plaquette loop. We have 



(4.1) 



[{T,UTi)pi,{T,Vr^)p^\ = 

(IV UTi)^i, TV VTj]pj - (TV yr,)[p,-, TV Un^ + (TV C/ri)(TV VTj)^i,Pj] 

and want to investigate what happens when the momenta are ordered to the left, i.e. whether 

[p,(IVC/ri),p,(TVFT,)] -p,(IVC/ri)[pi,TVFr,-]+ 
Pi{T^VTj)^j,T,Un] - %,P3\{T^Un){TrVTj) 

vanishes or yields new terms of order fi. Rewriting pj(Tr C/xj) = {Tr UTi)pi + [pj,TrC/rj] and 
using equation (4.2), expression (4.3) becomes 



[(TV Un)pi, \pj,Tr Vtj]] + [\pi, TV Un], (TV VTj)pj]- 

\pj,{TTUTi)[pi,TrVTj]] - \Pi,{Ti:VTj)[TrUTi,p,]] - (Tr C/ri)(Trl/r,-)] = 

(Tr Un)[p,, [p„Tt Vt,]] + (Tr yT,)[[p., Tr Ut,],p,] - (Tr Un)[p,, [p„ Tr Vr,]]- 

(TVyr,-)[Pi, [TvUn,pj]] - [\Pi,pj],TTUn]iTtVTj) - {TtUn)[\Pi,pj],TTVTj] = 

9 



(4.4) 



where we have made repeated use of the basic commutators (2.1). The last equahty in (4.4) 
holds by virtue of the Jacobi identity satisfied by the basic quantum operators. This latter 

property is not entirely trivial. For example, the holonomy and momentum operators as they 
arc usually defined in the continuum loop representation do not satisfy a Jacobi identity [22,9], 
a feature that gives rise to the occurrence of anomalous commutators for certain composite 
operators. 

It follows immediately that also for an arbitrary factor-ordering of the lattice constraints, 
a(Tr UTi)pi + (1 — a)pi{TrUTi), < a < 1, no anomalies appear. This is in particular true 
for the case a = ^, where the regularized constraint operators are self-adjoint. 



5 Summary 

We have performed a computation of the commutator of two regularized diffeomorphism 
constraint operators in lattice gravity, and found that their algebra closes without anomalies 
in the limit of vanishing lattice spacing. The discretization and quantization of the classical 
diffeomorphism phase space functions is not unique, but our result is independent of the 
choice of a local symmetrization and factor-ordering. Independence of factor-ordering followed 
from the simple structure of the constraints (linearity in momenta) and the fact that the 
basic lattice operators satisfy the Jacobi identity. We do not expect a similar behaviour for 
commutators involving also the Hamiltonian constraint. 

There are not many regularization schemes for full four-dimensional quantum gravity in 
which a computation of this type could be performed. Within our lattice formulation, it would 
in principle be preferable to have an exact remnant of the diffeomorphism symmetry at each 
stage of the discretization. This would enable one to study invariant measures and quantum 
states before the continuum limit is taken. One way of proceeding is to try to identify a dis- 
crete subgroup of the diffeomorphism group (rather than using the discretized "infinitesimal" 
generators) for each finite lattice, an issue we are currently considering. Transformations of 
this type will presumably be non-local in terms of the lattice variables. 
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Appendix 



In this appendix we present some intermediate results of the commutator calculation, 
depending on the symmetrization chosen. This will be helpful to anybody attempting a 
similar calculation. We split up the contributions into those that arise from contributions 
based at the same lattice point (the origin in our case) and those involving contributions 
based at neighbouring points on the lattice. It is understood that A, F and E are operators. 

(a) No symmetrization; contributions at the origin: 

~a^{ NH'fFi^iE] + E'i + Ef) + N'M^FI^ + F^^ + ^3^)^^ )- 

{M'^N^ + N^M'^)Fi2{E} + Ef + E^) + (M^AT^ + N^M^){Fl2 + Fjg + F^i)Ef+ 

^N^M^iVi + V2)Fi2{El + Ef + Ef)+ 

^N^M\{Vi + V2)Fi2 + (V2 + Vs)F^^ + (Vi + Vs)F^,)Ef+ 
N^M'^Fi^iViEl + V2Ef + VaEf) + N^M^iFl^ + F^^ + Fi^)V^El ) + O(a^) 

(Al) 

contributions involving neighbouring vertices: 



a^iN^M^Ff^iEj + Ef + Ef) + N^M^iFf^ + FJ3 + F.ii)Ef ) + 

{M'^N^ + N'^M'^)Fi2{El + Ef + Ef) + {M'^N^ + N'^M'^){F12 + F^^ + F^i)Ef+ 

^N^M\diFl2 + Se'jkA{F^2 + ^2Fi2)El + ^N^M\d2Fi2 + Se'jkAiF^^ + ^iFi2)Ef+ 

^N'M\2{Vi + V2)Fl2 + (V2 + V3)Fj3 + (Vi + V^)F^^)Ef- 

N^M^Fi^iViEf + V2Ef + VsEf) - N^M^Fl^idiE] + diEf + ^3^^ )+ 

N^M\Fl2 + Fi^ + Fi^)V^Ef - N^M^Fl^dsEf + F^^diEf + F^^d2Ef)- 

N^{diM^)Fi2Ej - M^{d2N^)Fi2Ef - N^{d^M'^)Fi2Ef - M^{dzN^)Fl2Ef- 

M'^{diN^)Fi^Ef - N\d2M^)F^iEf + N^M^eijki^AlF^^ + ^2^^ + A{F^^)E^' ) + 0{a^) 

{A.2) 

(b) Symmetrization over shift functions; contributions at the origin: 
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-a\^N'M^Fi^iEj + + 2Ef))- 

^(M^iV^ + M^iV^ + 2N^M'^)Fi2El + h^2M'^N^ + N^M'^ + N^M'^)Fi2E}+ 

-M2(iV^ - N^)F^^Ef + -N^{M^ - M'^)F^^Ef+ 
^N^M\{V2Fl2)El + {ViFl2)Ef + ((Vi + Vs)^!^^;^ )+ 
^Fl^iViEl + V2^f + 2V3^f ) ) + O(a^) 

contributions involving neighbouring vertices: 

a^^N^M^Fi^iE} + Ef + 2Ef))+ 

^(M^iV^ + M^iV^ + 2N^M'^)Fl2El + ^{2M'^N^ + iV^M^ + N^M'^)Fl2E'j+ 

^(SM^iV^ + M^TV^ + STV^M^ + N'^M'^)Fi2Ef+ 

^M'^{N^ - N^)F^^Ef + ^Ar^(M2 - M'^)F^iEf+ 

^N^M\{V2Fi2)E} + {ViFl^)Ef + ((Vi + Vs)^/^)^^,^ )+ 

N^Mh,,k{AiF^2E^' + + ^Fi*2(Vi^,' + Va^? + 2V3^f )- 

N^M^Fl^idiE] + a2£;2 + Ss^f ) - N^M\Fl^d^Ef + Fa^g^ii^f + Fg^Sai^f )- 

N^{diM'^)Fl^El - M'^{d2N^)Fl^Ef - N^{d^M'^)Fl2El - M^{d3N^)Fl2Ef- 

M^{diN^)F^sEf - N\d2M^)F^iEf + N^M^€ijk{2AiF^2 + A^F^i + A{F^3)E^' ) + O(a^) 

(AA) 

(c) Symmetrization over momenta; contributions at the origin: 

a^{-^N'M^Ef{Fl-, + F3\) ) + a^-^iM^N' + N'M^)EfiFi^ + F^,)- 
-iViM2E3(V2F2^3 + ViFg^) - -Fl,iViEl + V2£;| + 2V3F3^ ) + 0{a^) 



contributions involving neighbouring vertices: 
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a\^N'M\F^^ + F^,)Ef)+ 

^{M^N^ + N^M^){Fls + F^^)Ef + N^M^eijkiA{E^' + A{)F^2E^')+ 

^N'm\V2F^, + ViF^,)Ef + ^Fl^iViEj + V2E^ + 2V3i^f )- 

NH'I^Fl2{diEj + d2Ef + 93i?f ) - N'M\Fl2dsEf + Fs^a^ii^f + F^^d2Ef)- 

N^{diM^)Fi2E} - M^id2N^)Fi2Ef - N^{d3M^)Fi2Ef - M^{d3N^)Fi2Ef- 

M\diN')Fi^Ef - N\d2M^)F^,Ef + M\i^k{24Fi2 + A{fI^ + A{Fi^)E^' ) + ©(a^) 

{A.Q) 

(d) Symmetrization over both shift functions and momenta; contributions at the origin: 



-Fi^iViE} + V2Ef + 2VsEf) ) + O(g^) 



contributions involving neighbouring vertices: 



iViM2e,,fe(4Ei^ + £2)F^2E'' + \f{2{ViEI + V2EI + 2V3£;f )- 
N^M''Fl2{diE] + 92^f + 93£;f ) - N'M\Fl2d3Ef + Fj35i£;f + F3\92£;f )- 
N^{diM'^)Fi2El - M'^{d2N'^)Fi2Ef - N^{d3M'^)Fl2Ef - M'^{d3N^)Fl2Ef- 

M'^{diN^)F^sEf - N\d2M'^)F^iEf + N^M^eijk{2A{F^2 + A{F^i + 4f^^)E^' ) + O(a^) 

(^.8) 

The more one symmetrizes, the more cancellations occur already individually among 
terms based at the origin and among terms coming from the more non-local commutators. 
For example, terms of order still appear in cases (a), (b) and (c), but no longer when the 
totally symmetrized lattice operators are used. This is in line with the general rule that for 
lattice discretizations the convergence properties in the continuum limit are improved when 
one symmetrizes the lattice expressions over positive and negative lattice directions. Note also 
that terms containing partial derivatives of the shift functions M and N (that are crucial in 
obtaining the correct commutator) come from the commutators involving neighbouring sites. 
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